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Abstract 

In this paper we consider the problem of calculating the quantiles of a risky position, 
the dynamic of which is described as a continuous time regime-switching jump-diffusion, by 
using Fourier Transform methods. Furthermore, we study a classical option-based portfolio 
strategy which minimizes the Value-at-Risk of the hedged position and show the impact of 
jumps and switching regimes on the optimal strategy in a numerical example. However, the 
analysis of this hedging strategy, as well as the computational technique for its implemen- 
tation, is fairly general, i.e. it can be applied to any dynamical model for which Fourier 
transform methods are viable. 

Key words: regime switching jump-diffusion models. Value at Risk, risk management, 
Fourier transform methods. 
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1 Introduction 



Quantiles computation of profit and losses of a given financial position is a basic task for mea- 
suring and managing portfolio market risks. The dynamic of the driving risk factors, as well 
as the chosen risk measure are the basic ingredient for the analysis of any risk management 
strategy. In this paper the model we consider for the risky position is of the form S{t) = SQe^^^\ 
where X{t) is specified on a filtered probability space (il, {Tt}, T, V) as a jump-diffusion whose 
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parameters change over time, driven by a continuous time and stationary Markov Chain on the 
finite state space S = {1, 2, . . . , M}, representing the unobserved state of the world. In fact, em- 
pirical studies on the behavior of financial markets show the ability of regime-switching models 
to capture some peculiarities in the observed data, as firstly highlighted in the seminal paper 
by Hamilton [T3]. Since then there has been a growing effort in applying switching models to a 
wide class of financial and/or economic problems, such as time series analysis, portfolio theory, 
derivative pricing and risk management. On the other hand, the necessity of including jumps in 
the underlying models to provide better representation of their dynamical properties is widely 
recognized (see e.g. [9]). Regime-switching jump-diffusions turns out to be an appealing and 
flexible class of dynamic models. 

Among many different risk measures proposed in the literature, Value-at-Risk, although 
sharply criticized for the lack of sub-additivity and its inability to quantify the severity of an 
exposure to rare events, has been adopted as a benchmark in the financial industry and for 
regulatory purposes. It plays a central role in banking regulation and internal risk management, 
mainly due to its simplicity. We therefore take the VaR as a starting point of our analysis of 
risk management strategies. The computation of VaR in Regime-Switching models has been 
considered by several authors mainly in discrete-time setting (see e.g. [6] or more recently |14] . 
|23j). Here we consider this problem directly in the continuous time framework: as a matter 
of fact, the required computations can be very efficiently implemented with the help of Fourier 
Transform methods (see e.g. [2^). The use of this kind of technique for the analytical calculation 
of VaR has been considered in Duffie and Pan [12] who exploited the classical characterization of 
the distribution function in terms of the Fourier inversion of its characteristic function. The use 
of Generalized Fourier Transform and the FFT algorithm is more recent: see Le Courtois and 
Walter [17] who calculate the VaR for the Variance Gamma (VG) model and Kim et al. |15j . 
Scherer et al. [22], who consider the class of tempered stable and infinitely divisible distributions. 

As an application to risk management, we investigate the influence of jumps and switching 
regimes on the exposure to an underlying risky asset. More precisely we study a classical hedging 
policy based on options followed by an institutional manager whose aim is to minimize the VaR 
of a position. This type of analysis has been initiated by [1] a decade ago for a portfolio made 
by a risky asset following a log-normal random dynamic, and hence analytically solved in a 
Black-Scholes setting. More recently, it has been considered for a bond portfolio in [11], [3], [1]. 

By taking the VaR as the risk measure for potential losses L of a portfolio at a given level a 
(i.e. the value VaRa(L) such that Prob{L > VaRa(L)} = a), the strategy considered consists 
on minimizing the VaR of the option-hedged portfolio L^'^ with respect to the strike price and 
the quantity of the put option written on the risky asset, subject to a budget constraint: in 
other words, we hedge the risky position by buying a fraction € (0, 1] of a put option with 
maturity T and strike price K, but what K and h? The optimal hedging strategy is therefore 
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given by the solution of the fohowing program: 

< 

, hUP{K,T) = C, hG{0,l] 

IIq{K, T) being the time put price with strike K and maturity T. To implement the strategy, 
it is therefore needed i ) the calculation of the VaR for the risky asset and ii ) the corresponding 
value of a put option Hq {K,T). Under the classical no- arbitrage assumption, the price of 
the put option can be represented as the expected value of the payoff with respect to a risk- 
neutral probability (see e.g. [7j). Conversely, VaR is obtained under the objective or historical 
probability measure. Both steps can be efficiently faced with the Fourier transform technique. 

The paper is organized as follows: we firstly derive the optimality conditions for the VaR 
minimizing strategy (Section 2) and then we specify the Fourier Transform technique for calcu- 
lating quantiles and put/call option prices in a very general setting (Section 3). In Section 4 we 
introduce the regime-switching dynamic model, its generalized characteristic function and the 
main change-of-measure result for switching from the historical to the risk-neutral probability. 
Finally, in Section 5 some numerical experiments are reported to show the impact of jumps 
and regime-switching on the process quantiles and on the optimal hedging strategy. We use in 
particular a simple two-state model with gaussian jumps and quantify in such a case what is 
the effect of a wrong model choice. 

2 VaR and optimal risk management 

To measure the risk of a financial position the quantiles of its distribution function are commonly 
used. Given a confidence level a G (0, 1), the set of e-quantiles of the random variable Y is the 
interval [<laO^)-,1a(Xy\ where 

q-{Y) = inf{g G R|P(y < q) > a}, q+{Y) = inf{g G W\P{Y < q) > a]. 

For a random variable having continuous and strictly increasing distributions function Fy {y) , 
1a.(X) — 1a^) = IctO^) s-iid qa{Y) = Fy^(a), i.e. it solves the equation 

P{Y < qa{Y)) = a. 

Here we take the portfolio loss L to describe a financial position in a fixed time interval and, 
in order to simplify notations, we assume in this section that L has a continuous and strictly 
increasing distributions function. The Value-at-Risk at level a is defined as 

VaR„(L) = inf{g G M|P(L < q) > a). 
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Let St be the value of the risky asset, t € [0, T] and r be the risk-free rate, that without loss 
of generality we consider fixed in the period: we define the loss at time of such a position as 

L" = So - e-^^^T. 

Then 

VaR,(L") = 5o-e~^^^gi_„(ST). 

Let us now consider a classical hedging problem in which an institution has an exposure 
to a risky asset St and decide to hedge such an exposure in the interval [0,T] by buying a 
fraction h € (0, 1] of an European put option on the asset with maturity T and strike price K. 
Analogously to the situation considered by Ahn et al. (1999), we take as the hedged position 
the portfolio composed by the risky asset and the put option: the loss of the hedged portfolio 
at time is therefore 

L^^^ = So + hU^iK, T) - e-^^(ST + h{K - 5t)+) 

where IVf{K,T) is the price of the put option at time t. By defining the strictly increasing 
function 

g{u) =u-h{u- K)+ + hU^{K, T) 
where K = Sq — e~^'^K, it is immediately seen that 

therefore 

VaR,(L'^'^) = <7(VaR,(L")) = So - e~'"'' qi^aiSr) + hU(;{K,T) - e"^^/i(K - qi-aiSr))^ 

= VaR„(L^) + hU^{K,T) - e-'^hiK - gi_«(5T))+. 
Let us firstly notice that ii K < qi-aiSr), then 

YaRaiL'''^) > VaRc,(L") 
since IIq{K,T) > 0. Therefore the optimal hedging strategy is given by the following problem: 
' min;^,,, VaR„(L«) + hUP{K,T) - e-'-^h{K - (?i_„(St)) 
< hIlP{K,T)=C, (1) 
^ /i€ (0,1], ir>Qi_„(5T), 
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C being the budget constraint. Since h = C /lif {K,T), the optimahty first order condition for 
K is given by the following non-linear equation: 

Ii^{K,T) = {K-qi.^{ST))^U^{K,T). (2) 

Assuming that ([2]) has a solution K* > qi-aiSr) and the twice differentiability of the price 
functional we can prove that this is actually a minimum since 

by the convexity of the price functional w.r.t. the strike. Correspondingly, the optimal amount 
of the hedging put option is 

h* = T. (3) 

We now assume the following: 

Assumption 2.1. The price of the put option can be represented as the discounted expected 
value of the payoff at time T under a risk-neutral measure Q: 

Uf{K,T) = e-''^E^[{K - 5t)+]. 

Furthermore, let Fs{s) = Q{St < s) be the cumulative distribution function (cdf) of the random 
variable St under such a measure: hence 

{K - s)+dFs{s) = e-^^ K / dFs{s) - / sdFs{s) 

'CO \ J —CO J —CO J 

(kQ{St <K)- j sdFs{i 



and 

^Yif{K,T) = e-^^Q{ST<K). 

We can finally prove the following property: 
Proposition 2.1. If K* > qi^c^iSr), then VaR«(Lo*'-^*) < VaRa(L"). 
Proof. Since K* and h* are characterized through ([2|) and (l3|), we get 

VaR«(L^*'^*) = VaRa(L") + C - e-''^h*{K* - qi-aiSr)) 
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= VaR,(L") + ^p^^^^^^ iK* - qi-aiSr)) (^nf (i^*,r) - e^^ 
From Assumption 12.11 we have 

Anf (K,r) - e-^^ = e-^^{q{ST <k)-i)<o. 



Therefore 

□ 



VaR,(L[; ) = VaR„(L-) + -p^— - gi_„(ST))e-'-^(Q(5T < A'*) - 1) < VaR,(L"). 



Remark 2.1. Notice that the optimality condition ^ under Assumption \2.1\ becomes 

r-K 



KQ{St <K)- f sdFsis)) ={K- (/i_,(5T))e-^^Q(5r < K) 

J —oo / 



which simplifies to 

^^^^£.dF«W=,._„(5ri (4) 

and depends on both the objective and the risk neutral distributions V and Q. Furthermore, it 
easily seen that the l.h.s. is equal to the conditional expectation ¥P[St\St < K] which is an 
increasing function of K bounded by'E,'^[ST\- Therefore, the eg. ^ has a unique solution if and 
only if qi.^{ST) <^^[St]. 

The minimum VaR as a function of the budget C is therefore 

VaR* = VaR,(y) + C ^ n(i^*,r) ' 

which is a hnear function with negative slope. In the plane budget-risk the previous equation 
describe an efficient frontier giving for each level of the budget C the minimum VaR. 

Remark 2.2. It is easy to show that the problem of looking for the hedging strategy with the 
minimum cost C{h,K) = hIl{K,T) for a target level v of VaR, results in the same first order 
optimality condition for K and that we get the same linear efficient frontier 

c = rn(A-.) = („ - vaR.(i-)) „,^.,/_'^;'._^^^^^,^ . 

in the budget-risk plane. 
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3 The Fourier transform method 



Fourier transform methods are efficient techniques emerged in recent years as one of the main 
methodology for the evaluation of derivatives. In fact, the no-arbitrage price of an european 
style contingent claim can be represented as the (conditional) expectation of the derivative 
payoff under a proper risk-neutral measure (see e.g. [7]). These methods essentially consist 
on the representation of such an expectation as a "convolution" of two Generalized Fourier 
Transforms. Since the value of most derivatives depend on a trigger parameter, two main 
variants have been developed depending on which variable of the payoff is transformed into the 
Fourier space. Here we consider the technique introduced in [5] which consider the generalized 
Fourier transform with respect to the trigger parameter. 

More formally, let 11(5', K) be the payoff at maturity of the derivative: for example, 11(5', K) = 
{K — S)^ is the payoff of the put option. The no-arbitrage price is therefore given by 

no = e-^^E2[n(5r,i^)]. 

Due to the exponential structure of typical underlying dynamics of the form St = sq exp(Xt), 
it is convenient to represent the payoff with respect to the new variables Xt = log(5r) — log(so) 
and k = log{K), in such a way 11(5^,-?^') = n(so exp(XT), exp(A;)) = II{Xt + log(so), log(ii')). 

Therefore, let us denote with n(x, k) an arbitrary payoff function and with Ilx{z) its gener- 
alized Fourier transform (GFT) w.r.t. k, that is 

U^{z) = [ e''''U{x,k)dk, z G C; 
under proper regularity conditions (see e.g. [H]), Fourier inversion gives 

U{x,k) = — e-'^'=n,(z)dz, 
in some strip of C, from which 

BtHo = E[U{X{T) + log(so), k)] = [ n(x, k)Qidx) 

JR 

1 riu+oo -I riu+oo /• 

- / e-''''Il,iz)dzQ{dx) = — / e-'^'^ / U,{z)Qidx)dz 



1 

2^ 



^-i.fcjES[nx(T)+log(so)(^)]rf^- 



' II/— oo 

In order to implement our program, we need to evaluate 
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1. the VaR of the hedged position: this step require to solve w.r.t. v the equation V{St < 
v} = a; 

2. the value of a put option n^(A:, T) = e-''V[{e'' - 5oe^r)+]_ 
Let us consider the following "payoff' functions, 

^>i(x, A;) = (e'' - e'^)+, and <^2{x, k) = I{^<fc}, 

in such a way n^(K,T) = e-''^E^[<^i{XT + log(5o),A:)] and V{St < v} = V{Xt < k} = 
W,[^2iXT, k)], with k = log(u/S'o). Their GFT w.r.t. the trigger parameter k are 

^i(x,z) = ^, u>l and ^2(x, z) = -e''''' , u>0 

IZ — Z'^ z 

giving therefore the formulas 



1 pw+oo pilog(5o)(z-i)^,^('y ;\ 

U^{k,T) = e~^^^ e-'^^^ —^p^dz, u>l 

^uk q1—v / />+oo 



TT v7o i/^ — — 1/ + in(l — 2z^) 



and 



V{St<v} = -1 r "e--i°^(^'/5o)M£)rf^, ^>0 

<j)x{z) being the GFT of the process under the appropriate measure. If this is a regular 
functions in a properly defined strip of C, the transform method can be applied in both cases 
(see Lee (2004)). 

Since under the Assumption ()2.1|) the optimality condition is 

ii^{K,T) _ n^(K,r) 



{K-qi_^{ST)) 



the optimal hedging strategy is then implemented by running 
1. root search algorithm to find the value q* solution of 

TT V./n u-m J 
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2. root search algorithm to find the value K* solution of 



{K* - q*) 



5o»(/o+°°e- 



ml 




) 



Numerical quadrature must be used for integral evaluation. Alternatively the FFT algorithm 
can be used to efficiently approximate integrals (see [18]) and then a standard root-finding routine 
will find the required solutions. 

4 Regime- Switching Jump Diffusions and measure change 

Let us consider on a filtered probability space {il, J^, TtjV) a stochastic process of the form 
St = 5'oe^% > 0, modeling the value, or price, of a risky asset for t e [0,r]. We consider a 
jump-diffusion setting in which the jump process is described as a marked point process (MPP), 
the parameters of which are driven by a finite state and continuous time Markov chain. 

We briefly recall here (see e.g. Runggaldier (2003)) that a MPP can be characterized through 
the couple (T^,!^), where {T^} is an univariate point process on M"*" and {Yn} is a sequence of 
random variables on a given measurable space {E,£), as a random measure p{dy,dt) for which 



Nt being the Poisson point process. The corresponding intensity for p is a measure-valued 
process \t{dy) for which 



is a martingale for each predictable process % and it characterizes the MPP. A common form 
of the intensity is Xt{dy) = Xtmt{dy), where At represents the intensity of the Poisson counting 
process and mt{dy) is a probability measure on the mark space {E,£) describing the jump 
component. Finally, the couple {Xt,'mt{dy)) is called the {V , J-'t)-local characteristic ofp{ds,dy). 
This setting has been introduced in the financial literature by Bjork et al. (1997). Although 
jump diffusion models can be described in somewhat different ways, the approach based on MPP 
turns out to be particularly useful for managing absolutely continuous change of measures. As 
a matter of fact, for our application we have to specify the dynamic model under both the 
objective (or historical) measure V and an equivalent risk-neutral (or pricing) measure Q. 

Let a{t) be a continuous time, homogeneous and stationary Markov Chain on the state 
space 5 = {1, 2, ... , M} with a generator H £ R^^^; furthermore, fi : S ^ R, a : S ^ R and 
^ : E X S ^ R are given functions, {E, £) being the measurable mark space. Without loss of 




n=l 
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generality, we can assume in the following C M. In a given interval < t < T, we consider 
therefore the dynamic 

dX{t) = {n{a{t)) - la^{a{t))dt + a{a{t))dW{t) + / 7(7/, a{t- ))p'' {dy , dt), X{0) = 0, (8) 
2 Je 

where W{t) is a standard brownian motion and p'^{dy, dt) is a MPP characterized by the intensity 

Xf{dy) = X{a)m{a, dy). 

Here A(-) represents the (regime-switching) intensity of the Poisson process Nt, while m(-, dy) are 
a set of probability measures on E, one for each state (regime) i € iS of the chain. The function 
7(2/, a) represents the jump amplitude relative to the mark y in regime a. Throughout the paper 
we assume that the processes a(-) and W{-) are independent and that W{-) and p°'{dy,dt) are 
conditionally independent given a{t). We denote = a{a{s) : < s < t} the cj-algebra 
generated by the Markov chain. Furthermore, we assume that E[e'''^-^'°^^] = e^^y''^^m{a, dy) is 
finite for each regime a. As usual, we also define the compensated point process q'^{dy,dt) = 
p°'{dy,dt) — X{a{t~))m{a{t~), dy)dt in such a way H{y,a{s~))q°'{dy,ds) is a martingale 

in t for each predictable process H satisfying appropriate integrability conditions. 

An application of the generalized Ito's Formula gives the corresponding jump-diffusion SDE 
for the asset price 

= /i(a(t))dt + a(a(t))dVF(t) + /(e^(^'"(*-)) - l)p"(dy,dt) (9) 



E 



Sit-) 

= [fiia{t)) + X{a{t))]K{a{t))))dt + a{a{t))dWit) + / (e^^^'^^*")) - l)g°(dy, dt), 
with 5(0) = So and K{a) = E[e^(^'") - 1]. 



Measure changes. An absolutely continuous transformation of measures in a jump-diffusion 
setting allows to change the intensities of the MPP and the Markov chain in addition to the 
translation of the Wiener process (see Runggaldier (2003)). In this context it results convenient 
to represent the underlying Markov chain itself as a MPP (see Landen (2001)) in such a way 

da{t) = I 6{z)v{dt,dz) (10) 

where v{dt,dz) is a marked point process with finite mark space (S,P(S)), S = {2; = {i,j) : i 7^ 
j, i,j € {1, 2, ... , N}} and 'P(S) = 2^ and compensator 

Xy{t,a{t-),dz) = J2hiAiait-)=i)<^{i,j){dz) (11) 
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e(j being the Dirac measure. In the previous formula the numbers hij are positive and 
such that J2f:^i,j=i hij = -ha, for i = 1, . . . , iV. 

Consequently, let '!9j(a) be a square integrable predictable processes, ht{a,y) a non-negative 
function such that 

ht{a,y)m{a,dy) = l,Vt G [0,r] 



and let ip{a) and ^{z) be strictly positive functions defined on S and S, respectively. We can 
define a new measure Q on the measurable space by setting 

dW^ = dWt - ^t{a{t))dt 

q'^'°'{dt,dy) =p{dt,dy) - ip{a{t))\{a{t))ht{a{t),y)m{a{t),dy)dt (12) 
\v{t, a, dz) = (^{z)\^{t, a, dz). 

Besides the translation of the Wiener process W^, we perform a change in the intensity of the 
MPP giving the compensated process q'^'°'{dt, dy) with (Q, J^()-local characteristic ('0(a(t))A(a(t)), 
ht{a{t),y)m{a{t),dy)) and a change of the intensity of the Markov chain which under Q has 
generator = {hf.} where 

N 

hf^ = ^i,j)h,j, h2 = - J2 ^ii^k)hik, i = i,...,N. 

By taking the Radon-Nikodym derivative 



Lf = exp 



10 JO 

+ / (1 - ^(a(s))A(a(s))a!s + / f \og{iP{a{s))ht{a{s),y))p{ds,dy) (13) 

JO JO JE 

+ / {1 - ^{z))X^{t,a{t-),dz)ds + \og{^ {z))v{ds,dz) 

Jo Jo Js 

and supposing that E'^[Lt] = 1 for all t, we have a probability measure Q on equivalent to V 
with dQ = LxdV, under which 

= [fiiait)) + aiait))Mait)) + V'(«(t))A(a(t))K« (a(t))] dt + a{ait))dW^ (t) 



+ (eyivMi-)) -l)qQ'»{dy,dt), (14) 
Je 

where ^^(a) = EQ[e'>(^'") - 1]. 
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In order to price derivatives under the model ([9]) we need to specify a risk-neutral or martin- 
gale measure, that is a measure under which the discounted price process e~^^St is a martingale. 
This is done by taking 

r - /i(a) - ■0(a)A(a)K;Q(a) 

a = (15) 

a{aj 

from which we finally get the risk-neutral dynamic for the underlying 

= rdt + a{a{t))dW^{t) + / (e^^^^'^^*")) - l)g'3'"(dy, dt) (16) 

= [r-i;{a{t))X{a{t))K^{a{t))]dt + a{a{t))dW^{t)+ [ (e^^^^'^^*")) - l)p^'°((iy, dt). 

Je 

Consequently, from (jlSp and ([9]), the market price of risk is (see [21j ) 

p{a) = /i(a) — X{a)K{a) — r = X{a)[K,{a) — ^lJ{a)^^'^ (a)] — a{a)6{a) 



X{a) 



a{a)0{a). (17) 



27{J/,") - l)(l - ^p{a)ht{a, y))m{a, dy) 
Correspondingly, for the process X{t) we have 

dX{t) = [r-\a^{a{t))]dt + a{a{t))dW^{t)+ f ^{y,a{r))q^^^{dy,dt) 



[r - ^cT2(a(t)) - V(a(t))A(a(t))K^(a(t))]dt + a{a{t))dW^ {t) 

^{y,a{r))p'^''^idy,dt). (18) 



The measure transformation defined by (jl2p through (jl3p preserves the probability structure 
of the stochastic process X{t) under both V and Q. It worth noting that we can specify infinitely 
many equivalent measures Q. In practice, the usual way to select one of the equivalent measures 
is to calibrate the model to a set of observed data. 



GFT for regime-switching jump-diffusions. The formulas ([6]) and ([7]) depend from the 
dynamic model only through its generalized Fourier transform (j)x{z). In this section we report 
the GTF for the regime-switching jump-diffusion model. Since we have to consider the process 
Xx under two different measure, we derive its characteristic function for the following general 
dynamic 

dX{t) = ^{a{t))dt + a{a{t))dW{t) + / 7(y, a{t~))p{dy, dt) 

Je 

where 



r 



-a^{a) — il){a)X{a)K'^ (a) under the measure Q 



^^'^^ ' fi{a) — ^a^{a) under the measure V ^^^^ 



12 



the MPP p{dy, dt) has intensity 
X{a,dy) = 

and the Markov chain has generator Q = {qij}ij=i,...,N where 



'4}{a)\{a)h{a,y)m{a^dy) under the measure Q 
A(a)m(a, dy) under the measure V 



e*(*'j)/i. • under the measure Q 



N 



e ^ riij unaer tne measure Sii ■ ■ ■ , ■ ^ /^.^n 

= I h,, under the measure V ^^J = ^^■■■^N^^^3,-^'^q^^ = -^<l^u. (20) 

In |20j (see also |10j ) it was proved the following 

Proposition 4.1. Let (f>j{z) = 'K[e^^"'^'^^^^'^^] be the generalized Fourier transform of the jump 
magnitude under the given measure. Then, by letting 

M^) = ^e(j) + l^z^cT^j) - iA(j) (</-,■ (z) - 1) (21) 

and '&i{z) = 'Qj{z) — 'dM{z), we have 

(22) 

where 1 = (1, . . . , 1)' e R^^^\ 1(0) = (I„(o)=i, ■ ■ • ,Ia{o)=Af)' e M^^''^ and Q' is the transpose of 
Q. 

Different models can be recovered with simple linear constraints on the full parameter set of 
our dynamics Q, This follows by noticing that if ^(z) = C^cr{i) = a, A(i) = A and 4>i{z) = 

<j){z) we are implicitly assuming a unique regime so recovering the well-known characteristic 
function of the (single-regime) jump-diffusion dynamic ^priz) = exp(z^ + ^iz^a^ — iX{(j){z) — 1)) 
which includes the standard geometrical Brownian motion (GBM) (A = 0) and the Merton 
jump-diffusion models (JDM). By letting Aj = in (I22p we get the regime-switching version of 
GBM (RSGBM) and finally the regime-switching jump diffusion model (RSJDM) with the full 
set of parameters 

^i,ai,Xi,hij, i,j = l,...,M. 

The evaluation of the characteristic function requires to compute matrix exponentials for 
which efficient numerical techniques are available (see Higham (2009)); conversely, the case 
M = 2 can be considered explicitly. The following can be proved (see [20] and the references 
therein) . 
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Proposition 4.2. Let yi^2 be the solutions of the quadratic equation + {qi + q2 — '^0)y — \9q2 = 
and 

qfW = ^ (e^^^(yi + ^1 + 92) - ei'2^(y2 + gi + (72)) 

q^(^) = ^ (e^^^(yi + 91 + 92 - i^) - ^^'^{yi + gi + 92 - i^)) . 

Then 

^M'^'] = l„(t)=iqr(0) + lait)=2^l{0) 

and therefore 

□ 

5 Numerical results 

In this section we report results about two kind of numerical experiments: 

1. calculation of quantiles for RSJD. The objective is to study the behavior of such a quan- 
tities by varying diffusions and jumps parameters. 

2. valuation of the optimal hedging strategy in the regime-switching jump-diffusion frame- 
work. 

We consider a two-state regime switching version of the jump-diffusion model with gaussian 
jumps. This is defined by choosing 7(2/, a) = y and two kinds of normal jumps, i.e. Y{i) ~ 
J\f{ai,bi) from which K{i) = E[(e^W - 1)] = e°'+^?/2 _ 1, ^ = 1^2. The two state Markov chain 

aft) € 5 = 11,2} has generator under the chosen measure Q = 1 ]■ Let ai,Xi > 

V 92 -92 y 

and fii,i = 1,2 be given parameters: the regime switching jump-diffusion Merton model is 
defined as ^ 

dX{t) = i{a{t))dt + a{a{t))dW{t)+ I [ yp''{dy,ds) (23) 

Jo Je 

where ^(a(t)) is given by p^ . X{t,a{t),dy) = X{a{t))(pa(t)iy)dy is the intensity process of the 
Poisson jump component and X{a{t)) € {Ai,A2}, a{a{t)) G {o"i,(T2}, /i(a(t)) € {/ii,/X2}, (piiv) 
being the density of a normal distribution A/'(aj, hi), i = 1,2. 

All numerical procedures were implemented in the MatLab© framework. 
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5.1 Quantiles and VaR Calculation 

Quantiles were computed by solving the equation 

-IT \Jo v-m ) 

with respect to q (see Section 3). A standard root-search algorithm was used together with the 
Gauss-Lobatto quadrature for approximating the integral. Few milliseconds were needed to get 
the required quantile on an Intel© Core 15. 

We calculate the quantiles of the RSJD models for different values of the parameters under 
the historical probability V . We start by considering the simple regime-switching geometrical 
brownian motion (Aj = = 6j = 0): as expected, the effect of switching drift and volatility 
results in a sort of mixing behavior between the corresponding GBM models, see figures ([1]), 
([2]), and such effect becomes more evident for growing time T . The frequency of Markov chain 
switching, driven by the generator Q, induces a quite different behavior, whose impact also 
depends on the value of T (figures (l3|) , dH . 

It is well-known ([9]) that introducing (gaussian) jumps in the GBM dynamic has a great 
impact on the tails of the distribution, see figure ([5]). As before, the Markov switching generates 
a mixing effect: quantiles curves of the RSJD models are between the curves of the corresponding 
JD models without regime switching. 

5.2 Optimal risk management 

In the second set of numerical experiments we face two main issues: (a) Risk Reduction: how 
much reduction of risk is obtained by implementing the optimal hedging strategy in the RSJD 
framework and in particular what is the impact of jumps and switching regimes on the main 
quantities driving the strategy? (b) Misspecified Modeling: what is the effect of a wrong model 
specification which discards regime switchings and jumps, when they are indeed present in the 
market, and consider the simpler GBM model? 

(a) Risk reduction and model impact. Inside each model - GBM, JD, RSGBM, RSJD - 
we show the behavior of the optimal hedging strategy obtained by changing the value of some 
relevant parameters: this corresponds to specify different levels of the market price of risk. In 
the figures ([9]), (fTOl) . (fTT| (fT2]) and (fT3l) the optimal strike K* and the corresponding value of 
the VaR are depicted together with the risk reduction percentage 

Va^L^j^ 
VaR«(L«) ■ 

It is apparent how the the strategy can be effective in reducing risk. On the other hand, 
it can be noticed that a perturbation of a single parameter results in a change of the optimal 
strategy. 
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Optimal strategy RSJD 


Optimal strategy GBM (a) 


Prsjd 


T=.5 
T=l 
T=3 


64.7442, 0.7197, 37.0189 
55.5928, 0.6165, 47.1767 
41.6851, 0.5294, 62.3356 


65.6191, 0.8433, 35.3880 (0.2905) 
57.1579, 0.7644, 44.4718 (0.2689) 
43.5664, 0.7382, 58.6379 (0.2254) 


0.0132 
0.0148 
0.0157 



Table 1: Optimal hedging strategy {K* , /i*, VaR*) under the simulated true model (first column), 
the fitted GBM model (second column - estimated volatility) and the corresponding value of 
I^RSJD- Here ui = .3, 02 = .05, Ai = 2, A2 = .8, ai = 0.0, 03 = 0.0, 61 = 0.08, 62 = 0.15, qi = 1, 
(72 = 0.2; furthermore r = 0.5%, a = 0.01 and the budget constraint is C = 0.1. 





Optimal strategy RSJD 


Optimal strategy GBM {a) 


Prsjd 


T=.5 
T=l 
T=3 


61.1841, 0.0581, 45.2341 
45.8347, 0.0833, 60.5069 
18.8056, 0.4015, 83.7630 


63.3076, 0.0816, 41.3746 (0.3137) 
52.7089, 0.0744, 52.6106 (0.3047) 
32.7103, 0.0797, 72.5986 (0.2930) 


0.0166 
0.0255 
0.0640 



Table 2: Optimal hedging strategy {K* , /i*, VaR*) under the simulated true model (first column), 
the fitted GBM model (second column - estimated volatility) and the corresponding value of 
Prsjd- Here ai = .3, 02 = .05, Ai = 2, A2 = .8, ai = 0.05, 02 = -0.3, hi = 0.08, 62 = 0.15, 
qi = 1, (72 = 0.2; furthermore r = 0.5%, a = 0.01 and the budget constraint is C = 0.01. 

(b) Misspecified Modeling. In order to explore the model sensitivity of the optimal hedging 
strategy, we implemented the following exercise. We firstly fixed a RSJD model by choosing 
a complete set of parameters. Then we generated a set of call/put prices on which we cali- 
brate the GBM model: hence we run the optimal hedging strategy obtaining Kq^j^j, ^hsM ^^"^ 
correspondingly the minimal VaR, VaR^^^.j. We finally calculated the probability 

Ursjd = P(L^'*''^' > VaR^BA,) 

under the RSJD model. Results are shown in Table ([T]) and ([2]). Notice that even when the 
optimal strategies are similar, the probability that the portfolio loss exceeds the (optimal) VaR 
is greater than the fixed level a = 0.01. Of course this behavior depends on the choice of the 
parameters, but the underestimation produced by a wrong model choice can be quite severe: 
see Table (l3|), where parameters from a real data set were used (see |20]). 

6 Conclusions 

In this paper we considered the problem of computing the quantiles of a risky position described 
by a regime-switching jump-diffusion dynamic model. The knowledge of generalized characteris- 
tic function for this class of processes allowed us to use the Fourier Transform methods to design 
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Optimal strategy RSJD 


Optimal strategy GBM (cj) 


Prsjd 


T=.5 
T=l 
T=1.5 


38.3721, 0.2497, 66.0564 
26.6034, 0.4103, 76.3270 
19.6884, 0.6506, 81.8069 


60.0168, 0.0785, 44.8655 (0.3479) 
49.4859, 0.0737, 55.8926 (0.3320) 
41.9632, 0.0741, 63.4963 (0.3291) 


0.1165 
0.1304 
0.1430 



Table 3: Optimal hedging strategy {K*, h* ,VaR*) under the simulated true model (first column), 
the fitted GBM model (second column - estimated volatility) and the corresponding value of 
I3rsjd- Here ai = .27, da = .13, Ai = 6.8, A2 = .8, ai = -0.13, 02 = -0.34, bi = 0.08, 
62 = 0.15, qi = 6.5, q2 = 0.002; furthermore r = 0.5%, a = 0.01 and the budget constraint is 
C = 0.01. 

an efficient algorithm for the calculation of quantiles. With this same technique, we analyzed 
a static hedging policy based on the constrained minimization of the VaR of the option-hedged 
portfolio. Numerical examples showed the impact of jumps and switching regimes on the optimal 
strategy in a two-regime, gaussian jumps framework and moreover the risk of a wrong model 
choice. 

Some final comments can be briefly outlined. Firstly, notice that different kind of jumps, as 
well as the number of regimes, can be readily considered in our computational framework, such 
as the double exponential Kou model ([E]). Furthermore, the analysis of the hedging strategy 
is fairly general, that is it can be applied to any dynamical model for which Fourier transform 
methods are viable, for example it can be extended to Variance- Gamma or Bates models. Finally, 
besides the choice of different dynamic models, it would be interesting to consider alternative 
risk measures, such as the Conditional Value at Risk (CVaR). This is certainly less commonly 
used in finance industry, but it is widely used in insurance industry being a coherent, convex 
and stable risk measure (see [5]). 
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Figure 1: Regime-switching GBM, starting state ao ^ I, qi = q2 = 5 and T = 1. Comparison with 
standard GBM. Analogous results for starting state ao = 2. 



o,=Og=0.05, a^=^ 

1.3 r- 
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Figure 2: Regime-switching GBM, starting state = 1, 9i = <Z2 = 5 and T = 1. 



20 





Figure 5: Jump diffusion model with gaussian jumps: red lines correspond to GBM. T = 1. 
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Figure 6: RSJD model with gaussian jumps: switching volatility (left picture) and drift (right picture). 
T = 1. 
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Figure 7: RSJD model: switching variance (left picture) and mean (right picture) of the gaussian jumps. 
T = 1. 




Figure 8: RSJD model: switching jump arrival times (left picture) and intensity of the Markov chain 
Qi (right picture). T — 1. 
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Figure 9: Risk reduction in GBM model: the continuous line corresponds to a model in which the risk 
is not priced, the dotted line to a variation of +25% of the drift and the dashed line to a variation of 
—25%. The other parameters are cr — 0.1, r = 0.05, T = 1 and 5*0 = 1. 




Figure 10: Risk reduction in JDM model: the continuous line corresponds to a model in which the risk 
is not priced, the dotted line to a variation of +25% of the drift - left figure) and mean jump (a - right 
figure) and the dashed line to a variation of —25%. The other parameters are cr = 0.1, r = 0.05, T = 1 
and 5*0 = 1. 
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Figure 11: Risk reduction in JDM model: the continuous line corresponds to a model in which the risk 
is not priced, the dotted line to a variation of +25% of the jump intensity (A - left figure) and jump std 
{b - right figure) and the dashed line to a variation of —25%. The other parameters are <t = 0.1, r = 0.05, 
A = 1, a = 0, T = 1 and 5o = 1. 
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Figure 12: Risk reduction in RS model: the continuous line corresponds to a model in which the risk is 
not priced, the dotted line to a variation of +25% of the drift (/ii,/i2 - left figure) and volatility ((Ji,cr2 
- right figure) and the dashed line to a variation of —25%. The parameters are cti = 0.15, CT2 = 0.05, 
r = 0.05, qi = q2 = I, T = 1 and = 1. No substantial changes are observed for this kind of variation 
of the parameters qi . 
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Figure 13: Risk reduction in RSJD model: the continuous line corresponds to a model in which the risk 
is not priced, the dotted line to a variation of +25% of the volatility (cri,cr2 - left figure) and jump std 
(&i, &2 - right figure) and the dashed line to a variation of —25%. The parameters are cti = 0.15, = 0.05, 
hi = 0.1, 1)2 = 0.2, r = 0.05, qi = q2 = I, T = 1 and S'o = 1. No substantial changes are observed for this 
kind of variation of the other parameters Ai = 1, A2 = 0.1, ai = 0, 02 = —.1, Qi = 1. 
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